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ABSTRACT 

The bipolar HH 212 outflow has been mapped in SiO using the extended configuration of the 
Plateau de Bure Interferometer (PdBI), revealing a highly collimated SiO jet closely associated 
with the H 2 jet component. We study at unprecedented resolution (0'.'34 across the jet axis) 
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the properties of the innermost SiO "microjet" within 1000 AU of this young Class source, to 
compare it with atomic microjets from more evolved sources and to constrain its origin. The 
SiO channel maps are used to investigate the microjet collimation and velocity structure. A 
large velocity gradient analysis is applied to SiO (2-1), (5-4) and (8-7) data from the PdBI 
and the Submillimeter Array to constrain the SiO opacity and abundance. The HH212 Class 
microjet shows striking similarities in collimation and energetic budget with atomic microjets 
from T Tauri sources. Furthermore, the SiO lines appear optically thick, unlike what is generally 
assumed. We infer Tkm — 50-500 K and an SiO/H 2 abundance > 4 x 10~ 8 — 6 x 10~ 5 for rvn 2 
— 10 7 — 10 5 cm~ 3 , i.e. 0.05 — 90% of the elemental silicon. This similar jet width, regardless 
of the presence of a dense envelope, definitely rules out jet collimation by external pressure, and 
favors a common MHD self-collimation (and possibly acceleration) process at all stages of star 
formation. We propose that the more abundant SiO in Class jets could mainly result from 
rapid (< 25 yrs) molecular synthesis at high jet densities. 

Subject headings: Stars: formation - Radio lines: ISM - ISM: jets and outflows - ISM: molecules - ISM: 
individual objects: HH212 



1. Introduction 

Millimeter interferometric studies of the L1448 
and HH211 Class sources have identified high- 
velocity SiO jet-like emission p ossibly related to 



the primary protostellar win d (Guillot eau et al. 
199a IChandler et al Il200ll Iffirano et al. 1 12006: 



Patau et al. ll2006tlGueth et al. 112007^ 7 A similar 
SiO jet has recently been discovered by Codella et 
al. (2007; hereafter Paper I) in the HH212 H 2 out- 
flow in Orion (d ~ 450 pc) using the new extended 
configuration of the PdBI. This study revealed 
highly collimated SiO emission with a close spa- 
tial and kinematic correspondence to near-IR H2 
knots, indicating that both are tracing the same 
molecular jet component. In addition, an inner 
pair of SiO knots with no near-IR H2 counterparts 
was identified at ±1'.'5 of the central source, with 
a radial velocity range pointing to a high degree 
of collimation. Continuum data at 1mm further 
reveal a compact, optically thick source probably 
tracing a disk of diameter < 120 AU. Similar con- 
clusions were reached by Lee et al. (2007) in a 
lower resolution, multi-species study with the Sub- 
millimeter Array (SMA). 

In this second paper, we further exploit the un- 
precedented resolution of 0'.'34 HPBW across the 
jet axis provided by the extended configuration of 
the PdBI to carry out the first comparison of the 
properties of the Class SiO "microjet" in HH212 
with atomic microjets from Class I /II sources ob- 
served on similar scales. We identify several simi- 
larities suggesting that the same collimation (and 
possibly acceleration) mechanism is at work in 



Class jets as in later stages. We also demon- 
strate that SiO is optically thick and close to LTE 
in the inner knots (as is not generally assumed) 
and discuss possible origins for the abundant SiO 
in Class jets, compared to their more evolved 
counterparts. 

2. Jet collimation and velocity structure 
2.1. Present SiO observations 

Figure 1 presents PdBI SiO 5-4 maps from Pa- 
per I of the inner jet knots, separated into three 
velocity intervals. The SiO microjet is extremely 
narrow, with a typical transverse FWHM ~ 0'.'4 
at all velocities, i.e. an intrinsic width of 0'.'2 = 
90 AU after correction for the PdBI HPBW of 
0'.'34 in the same direction^. 

Figure 1 also shows that, in both lobes, the re- 
gion of blue/red overlap is not coincident with the 
region of highest radial velocities, but is slightly 
trailing behind it by about 0'.'4. The lack of 
blue/red contamination towards the fastest gas re- 
quires that its motions are highly forward-directed 
with a semi-opening angle < 4° (see Paper I). The 
blue/red overlap at low velocities ~ 3-4 km s _1 
traces less collimated, slower material in the wake 
of the fastest gas. 



1 A slightly larger width of 0'.'35 was quoted in Paper I, where 
cleaning had not been optimized for the innermost jet re- 
gions. 
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2.2. Comparison with jets from more 
evolved sources 

The width of atomic jets from T Tauri Class 
II sources spans a relatively broad range, depend- 
ing on the brightness of bowshock wings driven by 
internal working surfaces. In Fig. 2, the intrinsic 
FWHM of the HH212 SiO microjet is compared to 
the broadest (DG Tau) and narrowest (RW Aur) 
atomic microjets from Class II sources studied so 
far using ground-based adaptive optics or HST 
(jDougados et al. Il200d fWoitas et al. 1120021 ). We 
find that the HH212 SiO microjet falls exactly in 
the same range as Class II jets on scales 500-1000 
AU. Similar results are found for the SiO jet from 
the HH211 Class source (width of 95-125 AU at 
distances of 300-600 AU; Gueth et al. 2007). 

Also indicated in Fig. 2 is the width of the 
HH212 jet at 50 AU from the source, ~ 40 mas 
= 18 AU, as inferred from th e bow shape of 



H2O maser spots within 100 mas (jClaussen et al 
19981 ). Again it is undistinguishable from that of 



atomic microjets at the same distance. We thus 
find no evidence of a higher jet collimation in Class 
sources compared to the T Tauri stage where 
only a thin disk is present, although the dense in- 
falling envelopes characterizing the Class stage 
would be capable of strongly reconfining a radially 
expanding wind (jDelamarter et al. 1120001 ). This 
definitely rules out collimation by external pres- 
sure gradients and requires that jets from young 
stellar objects are self-collimated by internal mag- 
netic stresses. The jet MHD collimation process 
appears to be the same at all phases, with all fast 
material confined within a beam diameter of about 
15-20 AU over a distance ~ 50 AU. 

We further note that the HH212 Class mi- 
crojet follows interesting scalings compared with 
Class II microjets concerning its energetics and 
kinematics, (i) The mass ejection to accretion rate 
in HH212 estimated by Lee et al. (2007) from CO 
emission farther out along the jet is 15% (scaling 
with Vjet/100 km s _1 ). This is similar to the ratio 
of 1 0% found for spatial ly resolved Class II jets 
(e.g. IWoitas et ai~l l2002h . (ii) The HH212 knot 



speed of 100-15 km s 1 is typically half that in T 



Tauri jets (e.g. iDougados et al. 1 l2000h . for a four 



times lower stellar mass of O.15M (|Lee et al 
2006). Hence the jet speed appears reduced in 



the same proportion as the escape speed from the 



central object. Such scalings would be consistent 
with the jet acceleration mechanism and launch- 
ing zone also possibly being the same at all phases. 
However, similar data in a larger sample of Class 
jets would be needed to confirm this conjecture. 

3. SiO abundance in the HH 212 microjet 

3.1. SiO line ratios and brightness temper- 
atures 

In order to constrain the physical conditions as- 
sociated with the SiO emission in the inner jet, 
we compared the J=2-l and 5-4 line intensities 
from Paper I. For proper comparison, the SiO(5-4) 
map, originally obtained with a 0'.'78 x 0'.'34 resolu- 
tion, was reconstructed at the lower resolution of 
the SiO(2-l) map (l'/89 x 0'.'94). Figure 3 plots on 
a main beam ( "MB" ) scale the reconstructed 5-4 
line profiles at the peaks of the inner SiO knots, 
and the ratio Tmb (5-4)/Tmb (2-1) as a function 
of velocity (bottom panels). It can be seen that 
the ratio is ~ 0.75-1.1 across the blue knot profile, 
and ~ 0.5-0.85 across the red knot profile. Rela- 
tive calibration uncertainties between the 2-1 and 
5-4 lines are estimated to be ~ 20%. We similarly 
evaluate the SiO (8-7) to (5-4) intensity ratio by 
degrading our PdBI map to the 0'.'96 x 0"69 SMA 
beam of Lee et al. (2007). The resulting (5-4) 
spectra towards the inner SiO knots are also plot- 
ted in Fig. 3. Comparison with Fig. 10 of Lee et 
al. (2007) yields an (8-7)/(5-4) ratio in the range 
0.7-1 at all velocities. The relative calibration un- 
certainty could reach 30%. 

A third constraint is provided by the peak main 
beam temperatures Tmb (5-4) — 25 K in both 
knots in our original PdBI beam (top curve in Fig- 
ure 3). As the jet is broadened by a factor ~ 2 
by beam convolution across the jet (cf. Sect. 2), 
the intrinsic line temperature Ta(5-4) is at least 
25x2 = 50 K. Including beam dilution along the 
jet axis with O'.^S HPBW, the intrinsic line bright- 
ness could reach 200 K if the knot is roughly cir- 
cular. 

3.2. LVG modelling: evidence for optically 
thick SiO 

The line ratios and (5-4) intrinsic brightness 
are compared with the result of a large velocity 
gradient (LVG) code, which considers the first 20 
levels of SiO and the rate coefficients for collisions 
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with H 2 reported by Turner et al. (1992) up to 
T k = 300 K. We explored H 2 densities from 10 5 to 
10 7 cm -3 (see Sect. 3.3) and an LVG optical depth 
parameter n(SiO)/ (dV/dz) — Nsio/AV ranging 
from 10 12 to 10 17 cm -2 (km s -1 )" 1 , i.e. from the 
fully optically thin to optically thick regime. Our 
typical model results are illustrated graphically for 
Tk = 100 K in Fig. 4, and compared with observed 
values in HH212. 

We find that the usual approach of assuming 
optically thin emis sion to deriye 11CR2) and Tk 
from line ratios (e.g. Gibb et al. 112004 : iNisini et al. 
2007h would give inaccurate results in our case: 
As shown in Fig. 4(bottom panel), no optically 
thin model (starred symbols in the curves) can 
simultaneously reproduce the observed values of 
both SiO(8-7)/(5-4) and SiO(5-4)/(2-l). Values 
~ 1 for both ratios are only achieved when ap- 
proaching the optically thick LTE regime (Tr ~ 
Tk), which is the point of convergence of all den- 
sity curves at sufficiently high opacity. Wc infer 
that Nsio/AV must be greater than ~ 10 15 cm~ 2 
(km s" 1 ) -1 , while n(H2) is not well-constrained. 
The high (5-4) intrinsic brightness of 50 K-200 K 
also independently argues for a large optical depth 
parameter (Fig. 4, top). It also indicates that Tk 
lies in the range 50-500 K, or else the predicted 
7r(5-4) close to LTE would be too low/high. 

We note that substantial SiO optical depth 
could be rather common in the innermost part of 
Class jets, if they are as narrow as in HH212. 
In the L1448 jet, for example, a column density 
~10 14 cm~ 2 has been derived from single-dish 
measurements of the 5-4 transitio n assuming a 



ten: 



X(SiO) = 4x10" 



2007h . A narrower 



jet width of 2" (jNisini et al. 
width of ~0.2" would result in a column density 
higher by an order of magnitude, implying, as 
in HH212, a line optical depth larger than unity. 
Hence the low Tmb(5-4) ~ 0.1-1 K in single-dish 
observations could result mainly from severe beam 
dilution of the SiO emission, as argued previously 
by Gibb et al. (2004) and illustrated in Fig. 3. 
The SiO abundances would then be substantially 
larger than previously reported. 

3.3. SiO abundance and H2 density 

Noting that N Si0 /AV = n(SiO)/ (dV/dz), the 
SiO abundance with respect to H2 may be writ- 



10 15 cm- 2 km- 1 s 



10" 



n(H 2 



dV/dz 



410- 11 f 
(1) 

The adopted line-of-sight velocity gradient dV/dz is 
typical of cooled regions with Tk < 100 K at the rear 
of planar C-shocks and is probably a lower limit. A 
steeper gradient dV/dz is given by the ratio of the 
FWZI of the SiO line profile (~ 10 kms" 1 ) to the knot 
width (0'.'2 = 100 AU), which would increase X(SiO) 
by a factor 16 from the above formula. 

The main uncertainty in X(SiO) stems from the 
unknown H2 density in the SiO knots. A reasonable 
range may be inferred from the presence of shock- 
excited H2O masers at 0.1" from the source. Mag- 
netic field strengths and line ratios in H2O masers 
around YSOs typically require p reshock H nuclei den- 
sities n H ~ 10 7 - 10 s cm" 3 (|Kaufman fc Neufeld I 
1996). Assuming that density roughly drops with dis- 
tance as 1/r 2 (cf the DG Tau jet; Lavalley-Fouquet et 
al. 2000), one infers a preshock density ~ 10 5 — 10 6 
cm -3 at the SiO knots. Shock compression could in- 
creas e these values by a about a n order or magnitude 
(e.g. iKaufman fc Neufeld 1 1 19961 ) . so that the density 
is in the range ~ 10 5 — 10 7 cm -3 . The resulting min- 
imum SiO abundance for optically thick emission is 
X(SiO) > (4 x 10" 8 - 4 x 10" 6 ) x (1 - 16), with the 
higher value corresponding to the lower density, and 
the additional factor 1-16 arising from the uncertainty 
in velocity gradient. Assuming a solar abu ndance of 
(Si/H) ~ 3.5 x 10 -5 (|Grevesse et al. Ill998l h between 
0.05% and 90% of the elemental silicon is in the form 
of SiO. 

3.4. Origin of the SiO component 

Our PdBI observations of the HH212 microjet set 
stronger constraints than previously on the origin of 
the SiO in protostellar outflows, because of the shorter 
timescales involved and the unusually high collimation 
and SiO column densities indicated by our data. 

Given the proper motions of 60-150 km s -1 for 
H2O masers and H2 kno ts l|Claussen et al. 1998: 
iMcCaughrean et al. hooj ), the dynamical time of in- 
ner SiO peaks at 500 AU is only 25 yrs. SiO should 
thus be incorporated very rapidly in the flow. The for- 
mation of SiO in outflows is usually attributed to sput- 
tering of Si atoms from charged grains in a magnetized 
Gshock with ion-ne utral drift speeds > 25 km s _1 
(|Schilke et al. Ill997l ). Updated C-shock models with 
improved sputtering yields, SiO formation rates, and 
molecular cooling (Gusdorf et al., in preparation) show 
that the required conditions for optically thick emis- 
sion are reached for shock speeds 35-45 km s^" 1 and 
preshock densities of 10 5 — 10 6 cm -3 but only at the 
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rear of the shock where velocity gradients are small, 
i.e. after 400-150 yrs. As this exceeds the knot dy- 
namical time, non-steady truncated C-shocks need to 
be considered to model SiO-emitting shocks on such 
small scales. 

Another long-standing issue is whether SiO molecules 
originate from shocked ambient material or trace the 
primary jet itself. The SiO microjet diameter of 100 
AU is comparable to the centrifugal disk diameter of 
120 AU indicated by our 1mm co ntinuum size (Pa - 
per I) and by envelope kinematics l|Lee et al. Il2006h . 
Hence we would expect little infalling molecular mate- 
rial left on-axis to refill the jet path between successive 
ejection episodes, unless this material is very warm. 
This would appear to favor an origin of the SiO in 
the jet itself. The option is appealing, as the higher 
densities of Class jets, and the accompanying low 
temperature and ionization, are indeed conducive to 
molecular formation. In an early study of chemistry in 
protostellar winds, Glassgold et al. (1991) found that 
Si atoms are quickly converted into SiO at high mass- 
loss rates Mj e t > 10~ 6 M© yr" 1 . For a dust-free wind, 
the predicted SiO abundances are ~ 50-100% of the 
total elemental silicon. However, the recent finding of 
a substantial deple tion of Fe and Ca at the base of 
several Class I jets l|Podio et al. Il2006h indicates that 
jets are not dust-free and that grains are only partly 
eroded along the flow. In the HH34 jet, 13% of Fe 
has been returned to the gas at distances > 1500 AU. 
The same process at work in Class jets would release 
Si atoms in a sufficient amount to produce optically 
thick SiO emission if M jot > 1O~ 6 M yr~\ In the 
inner SiO knots of HH212, this mass-flux is achieved 
for n(H 2 ) > 10 6 cm~ 3 (with Vj ct = 100 km s" 1 and a 
jet radius of 50 AU), thus only 0.5%-9% of Si would 
be needed, if all is converted into SiO (see Eq. 1). 

4. Conclusions 

Our finding that jet collimation in the HH212 Class 
source is similar to that in T Tauri stars favors 
a collimation mechanism independent of the pres- 
ence of a dense envelope, i.e. most probably inter- 
nal MHD stresses. The ejection/accretion ratio and 
the jet speed/escape speed ratio also appear to be 
similar to those in Class II, possibly suggesting the 
same acceleration mechanism as well. The main dif- 
ference between Class jets and their more evolved 
analogs would then be their differing chemical compo- 
sition, with abundant molecules at the Class stage, 
a mix ed atomic-molecular com position at the Class I 
stage (|Davis et al. Il200ll . l2003h . and a purely atomic 
flow at the Class II stage. 

We also find that SiO is optically thick, so that its 
abundance is larger than previously estimated. The 



extremely narrow width of the SiO jet revealed by 
PdBI further argues that this species is not formed 
in swept-up material, but more likely within the jet 
itself. We thus propose that the higher SiO content 
of Class jets could mainly reflect an increase in jet 
density (hence, a higher efficiency of molecular forma- 
tion), linked to the increased mass-accretion rate at 
earlier stages. 
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0.5<|V-V sys |<5.5 km/s 5.5<|V-V sys |< 10.5 km/s |V-V gys |> 10.5 km/s 




10-1 10-1 10-1 

Offset (") 

Fig. 1. — SiO (5-4) emission maps of the HH212 microjet in 3 different velocity ranges. Blue and red 
contours refer to blueshifted and redshifted gas, respectively. A cross marks the position of the continuum 
source from Paper I: a(2000) = 05 /l 43 m 51 s .41, (5(2000) = -01° 02' 53"160. Contour spacing is 50 mJy/beam 
km s _1 with the first contour at 100 mJy/beam km s^ 1 . 
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Fig. 2. — HH212 intrinsic jet width compared to the range spanned by atomic microjets from Class II sources, 
corrected for the instrumental PSF (small dots: Woitas et al. 2002; large dots and triangles: Dougados et 
al. 2000); Our SiO PdBI measurements are shown as filled squares; the H2O maser width from Claussen et 
al. (1998) as a filled star. 



8 



1 

0.5 F 




HH212 Si0(5-4): North 2 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




i .... i . 



-L 



-15-10 -5 5 10 15 20 
V LSR (km s" 1 ) 



30 
25 
20 



HH212 Si0(5-4): South 2 



_l I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 

! 0.78" x 0.34" (21°) n 

j- 0.96" x 0.69" (38°) I 
! - - - 1 .89" x 0.94" (22°) 

f 

"■ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 III 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 

-i 

i.. 

p J 


1 

; i : 

■■ i . . . ■ i ■ ■ ■ ■ i ■ . ■ . i i . . . ■ i ■ ■ ■ ■ i ■ ■ ■ . i ■■ 



-15-10 -5 5 10 15 20 
V LS r ( km s" 1 ) 



Fig. 3. — Top panels: Line profiles in SiO J=5 4 towards the inner SiO knots at various resolutions: the 
original PdBI beam (solid histogram), the SMA J=8-7 beam (dotted histogram), and the PdBI SiO J=2 
1 beam (dashed histogram). Beam PAs are listed between parentheses. Note the dramatic decrease in 
brightness temperature with increasing beam dilution. The vertical dashed line marks the ambient LSR 
velocity (+1.6 km s _1 ; Wiseman et al. 2001). Bottom panels: Line temperature ratio Tmb(5-4)/Tmb(2-1) 
at the resolution of the PdBI SiO J— 2-1 map, as a function of velocity. 
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Fig. 4. — Top: SiO line temperature ratio T(5— 4)/T(2— 1) versus intrinsic line temperature Tr(5-4) for LVG 
slab models at Tk = 100 K. Each curve corresponds to the labelled log(n(fi2)), with dots marking values of 
iVgio/AF increasing (left to right) from 10 12 to 10 15 cm -2 (km s^ 1 )^ 1 by factors of 10. Symbols with error 
bars illustrate the range in line ratio in the inner SiO knots of HH212 (including calibration uncertainties) 
and the range in Tj>(5-4) after correction for beam dilution. Bottom: Same as above for the T(5-4)/T(2-l) 
ratio versus T(8-7)/T(5-4). 



